The Extraordinary 'Superthin' Spiral Galaxy UGC 7321. I. Disk Color 
Gradients and Global Properties from Multiwavelength Observations 



L. D. Matthews^'^'^ 
J. S. Gallagher, III^-^.s 
W. van DrieP'^'^ 

ABSTRACT 

We present B- and i?-band optical imaging and photometry, Ha narrow-band 
imaging, near-infrared H-hand imaging, and Hi 21-cm spectroscopy of the nearby 
(V/j=407 km/s), Sd spiral galaxy UGC 7321. UGC 7321 exhibits a remarkably thin 
stellar disk with no discernible bulge component. The galaxy has a very diffuse, low 
surface brightness disk, which appears to suffer relatively little internal extinction in 
spite of its nearly edge-on geometry (i ~88°). If seen face-on, UGC 7321 would have 
an observed central i?-band surface brightness of only ~23.4 mag arcsec"^. 

The UGC 7321 disk shows significant B — R color gradients in both the radial 
and vertical directions: A(i? — R) >0.80 magnitudes along the galaxy major axis, 
and A(B — R) as large as 0.45 magnitudes is observed parallel to the galaxy minor 
axis. These color gradients cannot be explained solely by dust and are indicative of 
changes in the mix of stellar ages and/or metallicity as a function of both radius and 
height above the galaxy plane. The outer regions of the UGC 7321 disk are too blue 
to be explained by low metallicity alone {B — R < 0.6), and must be relatively young. 
However, the galaxy also contains stellar populations with B — R >1.\, indicating it is 
not a young or recently-formed galaxy. 

The disk of UGC 7321 is not a simple exponential, but exhibits a light excess at small 
radii, as well as distinct surface brightness zones. Despite its organized disk structure, 
many of the global properties of UGC 7321 {Mb = -17.0; A4i^/=l.lx lO^A^©; 
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j^"' =1.1 Mq/Cq] W20= 233 km s~^; hj.=2.1 kpc) are reminiscent of a dwarf galaxy. 
Together the properties of UGC 7321 imply that it is an under-evolved galaxy in both 
a dynamical and in a star-formation sense. 

Subject headings: galaxies: spiral — galaxies: fundamental parameters — galaxies: 
evolution — galaxies: individual (UGC 7321) 



1. Background 

Ogorodnikov (1958) and Vorontsov-Vel'yaminov (1967) were among the first to draw attention 
to a subset of edge-on spiral galaxies with disk axial ratios as extreme as 30:10 and having little or 
no discernable bulge component. Goad & Roberts (1981) dubbed these galaxies "superthins" and 
recognized that spiral galaxies selected on the basis of their thin stellar disks and extreme axial 
ratios also tend to possess a variety of other intriguing properties. 

For example, Goad & Roberts (1981) discovered that superthins have low emission-line 
intensity ratios ([Nil] /Ha and [N ii]/[S ii]), which in a modern perspective suggests that superthins 
are moderate ionization, low metallicity galaxies more like irregulars than spirals in terms of their 
Hii region properties (see also Bergvall & Ronnback 1995). Similar trends were confirmed for 
larger samples by Karachentsev & Xu (1991) and Karachentsev (1991). 

In spite of having "classic" double-horned global Hi profiles (e.g., Matthews &: van Driel 
1999), often superthins exhibit slowly-rising, solid-body rotation curves throughout most or all of 
their stellar disks (Goad & Roberts 1981; Karachentsev & Xu 1991; Karachentsev 1991; Cox et al. 
1996; Makarov et al. 1997,1998; Abe et al. 1999). Such rotation curves are generally characteristic 
of late-type dwarf and irregular galaxies rather than normal spirals (cf. Casertano & van Gorkom 
1991). 

Not only are superthins fascinating objects in their own right, but the very nature of 
these galaxies makes them extremely valuable objects for exploring a variety of fundamental 
astrophysical problems, ranging from disk formation and evolution scenarios, to the interpretation 
of high-redshift galaxy populations, to constraining the stability and dark matter contents of disks 
(e.g., Karachentsev 1989; Zasov et al. 1991; Dalcanton & Schectman 1996). In spite of this, to date 
only a handful of detailed studies of individual superthins have been undertaken (e.g., Gallagher 
& Hudson 1976; Bergvall & Ronnback 1995; Cox et al. 1996; Matthews 1998; Abe et al. 1999). 

Measurements of the vertical structure and disk color gradients are of particular interest in 
the case of superthin spirals, as these can provide important clues to their evolutionary history. 



These disk axial ratios were measured on blue photographic survey plates; however, we find the axial ratios tend 
to become less extreme when measured at a fixed isophote on deeper, red-sensitive CCD images. For this reason, 
measurements of vertical disk scale heights offer a more robust means of quantifying the thickness of edge-on galaxies. 
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Vertical structure measurements provide vital constraints on disk stability, dynamical heating 
mechanisms, and dark matter contents (Matthews 1998; Matthews 1999, hereafter Paper II), and 
when combined with measures of global properties and disk color gradients, information on the 
star-formation histories and dynamical evolutionary histories of the galaxies can be gleaned. 

UGC 7321 is a superthin galaxy ideally suited to such investigations. Since it is nearby 
(D~10 Mpc; Matthews et al. 1999b), its structure can be reasonably resolved from the ground. 
In addition, UGC 7321 appears to suffer little from internal extinction (see below), and it is 
nearly exactly edge-on {i ~88°), hence projection effects will have a minimal influence on derived 
vertical structural parameters (e.g., de Grijs, Peletier, &: van der Kruit 1997). Finally, the analysis 
of UGC 7321 is simplified compared with other edge-on spirals due to the absence of a bulge 
component. 

Here we present measurements of the global properties of UGC 7321 based on new 
multiwavelength observations including B- and i?-band imaging and photometry, narrow-band 
Ha imaging, near-infrared H-hand imaging, and Hi 21-cm pencil beam spectroscopy. We also 
examine Hi aperture synthesis data obtained by Rots, Roberts, & Goad (private communication). 
We use our datasets to measure the radial and vertical color gradients in the disk of UGC 7321 
(Sect. 3) and briefly discuss what the color gradients reveal about the evolutionary history of this 
galaxy. In a second paper (Paper II), we present measurements of the vertical scale height of the 
UGC 7321 disk. There we combine vertical structural measurements with results from the present 
work in order to further constrain the stability and dynamical history of UGC 7321. 



2. Imaging Data 

2.1. Near-Infrared Imaging 

We obtained near-infrared (NIR) H-hand imaging observations of UGC 7321 in May 1997 
using IRIM on the 2.1-m telescope at the Kitt Peak National Observatory]]. IRIM employs 
a 256x256 HgCdTe NICM0S3 detector with l'!09 pixels, yielding a field-of-view of A'.Q per 
side. The detector has a gain of 10.5 e^/ADU and a readnoise of ~37 e~ per pixel rms. The 
observations were taken in the presence of moderate cirrus, hence our iJ-band data cannot be 
reliably photometrically calibrated. 

We obtained sets of 9 co-added, 20-s exposures of UGC 7321 at each of 35 telescope pointings. 
To insure accurate flatfielding and sky subtraction, we shifted the galaxy to a different position 
on the detector before each series of exposures. In addition, interspersed throughout the target 
observations, we obtained several series of "blank" sky exposures at regions adjacent to the galaxy. 



^Kitt Peak National Observatory is operated by the Association of Universities for Research in Astronomy, Inc. 
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Linearization corrections to the data were made automatically during data acquisition. 

Data reduction was accomplished using standard IRAFQ tasks. To remove the sky background, 
a median was made of several groups of temporally consecutive object and blank sky exposures. 
From these, a median-averaged combination of 25 dark frames was subtracted to remove the 
structure present in the detector dark current. Finally, the dark-corrected sky frame was 
subtracted from all dark-corrected object frames. 

Flatfielding was accomplished by dividing by a night sky flat composed of the median of all 
dark-corrected sky and object exposures. Lastly, the flatfielded object exposures were aligned and 
co-added. Our resulting H-hand image of UGC 7321 is shown in Fig. 1. Note the striking thinness 
of the stellar disk and the complete absence of a bulge component even in the NIR. At the disk 
center, the vertical exponential scale height is only ~140 pc (Paper II)! 

Fig. 2 shows the //-band data in the form of a contour map. Because the disk is so thin, 
the contours have been stretched in the z-direction for clarity. Note that the contours are quite 
regular, even near the galaxy center, suggesting that there is not significant dust absorption in the 
i/-band. 

In Fig. 3 we plot a major axis profile extracted from the //-band data. From this profile, we 
see that the disk brightness of UGC 7321 falls off somewhat more steeply on the east side (the 
righthand side of Fig. 3) than on the west side. It is unlikely this can be attributed to internal 
extinction effects (which should be small at NIR wavelengths), or to flatfield errors (which are 
quite small in our //-band images), hence this suggests a slight lopsidedness in the stellar disk of 
UGC 7321. 

Arrows overplotted on Fig. 3 denote "zones" between which the slope of the brightness profile 
change fairly abruptly. These zones can also be seem from visual inspection of both our NIR and 
the optical images and are an indication that the disk surface brightness of UGC 7321 deviates 
somewhat from a smooth, exponential distribution (see also Sect. 2.2.3). This makes it difficult 
to characterize the //-band major axis profile in terms of a single exponential scale length. As 
shown by Matthews & Gallagher (1997), this phenomenon is common among extreme late-type 
spiral galaxies, although the effect is generally lost if the light profile is azimuthally averaged. 
Matthews & Gallagher (1997) suggested these surface brightness zones may be a signature of 
limited dynamical evolution in some of the smallest and faintest spirals (see also Sect. 2.3.3). 



2.2. Optical CCD Imaging and Photometry 



*IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Associated 
Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science Foundation. 
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2.2.1. Data Acquisition 

We obtained photometrically calibrated B- and -R-band CCD imaging observations of 
UGC 7321 in June 1997 using the 3.5-m WIYN0 telescope at Kitt Peak. In addition, we obtained 
narrow-band Ha and uncalibrated S-and i?-band observations of the galaxy with the WIYN 
telescope in April 1997. 

The direct imaging camera on the WIYN telescope employs a thinned STIS 2048x2048 CCD 
with 0''2 pixels, yielding a field-of-view of ~7' per side. The CCD gain was 2.8 e~ per ADU, and 
readnoise was ~8 e~ per pixel rms. Exposure times were 900s in B, 600s in i?, and 1500s in 
Ha. Seeing during both the April and June observations was ~0"6 FWHM. The WIYN telescope 
takes advantage of the best seeing on Kitt Peak through a combination of thermal controls, active 
optics, and dome ventilation systems. This, combined with the small pixel scale of the WIYN 
CCD imaging camera, makes WIYN an ideal facility for obtaining the data necessary for detailed 
studies of disk structure and color gradients of galaxies. 

Our CCD images were reduced using standard IRAF tasks. Individual frames were overscan 
corrected and bias subtracted in the usual manner. Flatfielding was accomplished by dividing by 
a mean of 5 dome flats taken in the appropriate filter. The WIYN CCD contains a bad column 
which was corrected in our images by interpolation from the two adjacent columns. 

Both the April and the June observations were obtained during dark time. Photometric 
calibration for the June run was accomplished by monitoring two standard star fields from 
Landolt (1992) at two different air masses and then applying standard color term corrections to 
convert from WIYN natural magnitude system to the Landolt system. The total formal dispersion 
in our WIYN magnitude zero points is ±0.03 magnitudes; the validity of our transformations 
was checked via comparisons with other recent WIYN photometric solutions. Errors in the 
colors are approximately ±0.02 magnitudes. Sky brightnesses during our observations were 
f'B,sky ~22.1 mag arcsec"^ and ^iR^sky ~19.8 mag arcsec"^. 

Our i2-band image of UGC 7321 is shown in Fig. 4 & 5. Note the very diffuse, low surface 
brightness (LSB) appearance of the disk in spite of its edge-on geometry. Several background 
galaxies are clearly visible through the disk (e.g., Fig. 6), implying internal extinction is low. 

2.2.2. Ha Emission 

Our continuum-subtracted Ha image of UGC 7321 is shown in Fig. 7. Due to cirrus, we were 
unable to obtain a flux calibration for this frame. However, qualitatively, this image shows that 
Ha emission is present out to radii of at least r = ±2^5 (where r is the distance from the galaxy 
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center measured along the disk major axis in the plane of the sky). 

We see in Fig. 7 that the bulk of the Ha emission in UGC 7321 is confined to a region near 
the disk midplane roughly 4" thick, although diffuse emission is clearly visible at z-heights of 
up to ±8" at galactocentric radii \r\ <1'0. We note however that there is a relative dearth of 
emission within a region offset 8" to the cast of the disk center (Fig. 8). This "gap" is 8" wide 
and contains faint plumes of emission extending to z > ±7" out of the galaxy plane. We also see 
marginal evidence of several other thin filaments which are more extended in z-height. Deeper 
narrow-band images are clearly needed to establish whether the central Ha feature could be a 
signature of blowout from the disk, and whether UGC 7321 may indeed contain a diffuse, ionized 
medium extended to large z-heights, analogous to that seen in more luminous edge-on spirals (cf. 
Dettmar 1995 and references therein). 

2.2.3. Surface Photometry 

We performed surface photometry on UGC 7321 using routines from the IRAF STSDAS 
analysis package. With the "ellipse" program, we fit a series of 9 concentric ellipses to the 
UGC 7321 images. Position angle and ellipticity of the ellipses were determined from the 
outermost galaxy isophotes and were kept fixed throughout the fitting. The galaxy center was 
chosen to be the position of peak brightness; this location was identical in the B and R (and 
H) frames. Foreground stars, cosmic rays, and background galaxies within the aperture were 
removed via background interpolation using the IRAF task "imedit" . Sky values were determined 
by measuring the sky counts in each of several rectangular apertures at various locations on the 
image (see Matthews &; Gallagher 1997). 

Because UGC 7321 is so close to edge-on, elliptical apertures are not ideal fits to the galaxy 
isophotes. Nonetheless, these apertures allow us to measure aperture magnitudes and colors for 
UGC 7321, and permit an estimate of the mean exponential scale length of its disk. Although the 
ellipticity of the fitted isophotes is an unreliable method of determining the inclination of near 
edge-on galaxies, our derived value [i 88°) is consistent with the inclination we estimate from 
the slight asymmetry of the dust features along the galaxy major axis. 

Our derived photometric parameters for UGC 7321 are given in Table 1. Errors for the 
aperture magnitude measurements were computed following Matthews & Gallagher (1997) and 
take into account sky, fiatfield and Poisson errors, as well as the scatter in the photometric 
solution. Large-scale fiatfield errors are the dominant source of uncertainty; the maximum 
amplitude of the fiatfield variations was ~1% of sky in both B and R. Our aperture magnitudes 
are in good agreement with previous photoelectric values reported by Tully (1988; mB=13.86) and 
de Vaucouleurs &; Longo (1988; mR=12.99). 

Fig. 9 shows our azimuthally averaged B-band radial surface brightness profile of UGC 7321. 
Although we noted in Sect. 2.1 that the disk of UGC 7321 is slightly lopsided and shows deviations 
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from a perfect exponential, an azimuthally-averaged brightness profile permits determination of a 
"mean" scale length which is useful for offering a global characterization of the size scale of the 
disk. 

Because UGC 7321 is viewed close to edge-on, one must take into account projection effects 
in analyzing the projected light profile and in deriving an exponential scale length. For a disk 
with an exponential radial brightness distribution viewed exactly edge-on [i = 90°), the projected 
radial brightness profile along the galaxy midplane {z=0) is expressed as 

L{r) = Lo{r/hr)Ki{r/hr) (1) 

where Ki is the modified Bessel function of first order (e.g., van der Kruit &; Searle 1981). At 
small radii {r/hr « 1), this expression can be approximated as 

L(r) « Lo[l + (rV2/i?)ln(r/2/i,)]. (2) 

Note this implies a slight flattening of the light profile will be observed at small radii compared 
with a simple, unprojected e~^^^'' function. At large radii {r/hr » !)> one can write 

L(r) « Lo{TTr/2hr)^/'^ex.p{-r/hr) (3) 

(see van der Kruit & Searle 1981). Because of the y/r term, this produces a slightly less steep 
light profile than a pure exponential function with the same scale length. As a result, a scale 
length derived from fitting a simple, unprojected exponential profile to an edge-on disk will be 
overestimated. At r = 2-3hr this effect is ~10% (e.g., van der Kruit & Searle 1981). 

To derive the scale length for UGC 7321, we have fitted a model projected exponential light 
profile to our data. Our model is overplotted on Fig. 9. For the purpose of deriving a scale length, 
the slight deviation of UGC 7321 from the edge-on {i = 90°) case is not significant. From our fit 

wc derive a scale length of /ir,B=44" it 2 (compared with /ir,_B=51" ± 5" we derive from simply 
fitting the function L(r) = Loe"''/'^'' to the disk at intermediate r). 

A comparison between the data and our model (Fig. 9) shows that in spite of lacking a bulge 
component, UGC 7321 exhibits a light excess at small radii compared with the prediction of an 
exponential disk. The U-band measured central surface brightness of UGC 7321 (before correction 
to a face-on value) is 21.6 mag arcsec"^, while that predicted from extrapolation of the exponential 
fit to small r is ~0.35 magnitudes fainter. Interestingly, this excess is even more pronounced in 
the i?-band and appears to correspond to a distinct region in the color map of UGC 7321 (see 
Sect. 3). This strengthens the suggestion that UGC 7321 appears to have a multi-component disk 
that is more complex than a simple, pure exponential (see also Sect. 2.1). In addition, we note 
that at r > 120" our observed light profile falls off faster than the projected exponential model, 
suggesting the stellar disk of UGC 7321 may be truncated (cf. Barteldrees & Dettmar 1994). 
However, this latter trend should be confirmed with deeper observations on a wide-field CCD 
where sky subtraction and fiatfielding can be accomplished with greater accuracy. 
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2.2.4- Discussion: The LSB Nature of UGC 7321 

For an optically thin disk which is exponential in both the r and z directions, and which is 
observed inclined at 90°, the face-on central surface brightness will be hz/hr times the observed 
edge-on value. Here hz is the exponential scale height (cf., van der Kruit k, Searle 1981). 
Transformations from observed surface brightness values to face-on values at other radii may be 
derived from Equations 1-3 above. Using a thin-disk approximation (sec the Appendix) we have 
derived the additional corrections to these values required for the case where a is disk observed at 
an inclination slightly less than 90°. After also taking into account internal extinction corrections 
(see Sect. 3.2.2), we then find that if projected to face-on, the disk of UGC 7321 would have an 
observed 5-band central surface brightness of ~23.4 mag arcsec"^, an extrapolated central surface 
brightness of ~23.8 mag arcscc"^, and a mean total disk surface brightness fis ~27.6 mag arcscc"^. 
Thus UGC 7321 is a very LSB galaxy, and much of its disk would likely be nearly invisible if 
viewed closer to face-on. As shown below, the internal extinction in UGC 7321 appears to be quite 
low, thus its LSB appearance must result from a rather low current star formation rate. 

Further evidence that the seemingly "anemic" nature of UGC 7321 is due to minimal current 
star formation comes from comparing its blue luminosity to its far-infrared luminosity. Using 
the IRAS 60/xm and 100/xm fluxes for UGC 7321 from the NED database and Sage (1993), 
respectively, and following the prescription of Rice et al. (1988), we derive Li?/ii=7.8x lO^^C©, 
a value nearly two orders of magnitude fainter than the mean for Scd/Sd spirals in the UGC 
catalogue (Roberts & Haynes 1994). Using our blue luminosity (corrected for internal extinction), 
we find that Lpin/Ls ~0.08 for UGC 7321. A comparison with Rice et al. (1988) reveals 
that such low Lfir/Lb ratios are commonly found only in two classes of galaxies: (1) old, red 
early-type spirals, dE's and SO's with few young stars; (2) gas-rich, extreme late-type spirals with 
diffuse, extremely LSB disks. Examples of the latter group of objects include the nearby galaxies 
NGC 4395 (Sd IV), NGC 45 (Sdm IV), and IC 2574 (Sm IV-V). All three of these galaxies have 
Lfir/Lb <0.10, >1, blue optical colors, and highly transparent disks (e.g., Sandage 1961; 
Matthews et al. 1999a). 



3. Color Maps 

Using our B- and i?-band WIYN data, we have produced a, B — R color map of UGC 7321 
(Fig. 10). A variety of structure is evident in this map, including significant B — R color gradients 
in both the vertical and radial directions. We discuss these features in detail below. 

The global B — R color of UGC 7321 is 0.99, as measured within the outermost observed 
isophote (25.2 mag arcsec"^ in B before correction for inclination), and after correction for 
Galactic foreground extinction. This B — R color is a typical value for a normal, late-type spiral 
(cf. Lauberts &; Valentijn 1988; dc Blok, van der Hulst, & Bothun 1995). This suggests that 
in spite of its edge-on geometry, the disk of UGC 7321 does not suffer severely from internal 
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reddening due to dust. We return to this issue below. 

3.1. The Nuclear Region of UGC 7321 

Near the center of the UGC 7321 disk, our color map reveals a small, very red (S — i? Ril.5) 
nuclear feature, only a few arcseconds across (Fig. 10). This feature is offset ~5" to the east of the 
disk center as determined from the brightness peak of the R- and i?-band images. One possibility 
is that we are seeing the signature of an embedded nucleus. However, although compact star 
cluster nuclei are common in late-type low-luminosity spirals with diffuse, LSB disks (Matthews 
Sz Gallagher 1997; Matthews et al. 1999a), we do not see any direct evidence of such a feature 
in UGC 7321 in either our WIYN images or in images obtained with the Hubble Space Telescope 
(Matthews et al. 1999b). Nonetheless, the location of the compact red feature in the B — R color 
map is near to that of the peculiar Ha emission features shown in Fig. 8. The optical longslit 
spectrum of UGC 7321 obtained by Goad & Roberts (1981) also shows a possible kinematic 
disturbance in the Ha emission at this location. These are hints that some interesting physical 
processes are at work near the center of the UGC 7321 disk. Further investigation of this region 
via high resolution spectroscopy, deep Ha imaging, and NIR colors may be fruitful. 

Surrounding the compact red feature, our color maps reveal a more extended red region 
{B — R ~1.2), visible to r ±20" on either side of the disk center and showing a rather abrupt 
boundary (see also below). Intriguingly, the extent of this region corresponds very closely to 
the region over which we observe a light excess over a pure exponential disk fit (Sect. 2.2.3). 
This raises the possibility that this red central region might possibly represent an ancient central 
starburst, the core of the original protogalaxy, or perhaps even a kinematically distinct disk 
subsystem analogous to the bulge of normal spirals. 

3.2. Radial Color Gradients 
3.2.1. General Trends 

Cutting into the red central region discussed above are thin blue bands of stars visible along 
the midplanc of the galaxy. These bands grow both thicker and bluer with increasing distance 
from the galaxy center. At |r| =20", this layer has B — R Ril.05, reaching B — R RiO.85 at |r|=l!0, 
B-R ^0.80 at \r\ =1^5, B -R fa0.45±0.10 at r=2(7, and B - R ?«0.55±0.10 at r = -2(7. 
Thus, from the nuclear region (where B — R k, 1.5) to the outer disk edge we see a total B — R 
color change of up to 1.05 magnitudes along the major axis of UGC 7321. This is illustrated 
in Fig. llo, — d, which shows B — R color profiles of UGC 7321 extracted along the major axis, 
as well as 1"6 north (through the red nuclear feature), 1"2 south, and 4"6 south relative to the 
major axis, respectively. The extracted profiles were averaged over 12-pixel- wide strips and then 



-10- 



smoothed, hence the red color of the compact red central disk feature is slightly subdued in the 

profiles shown in Fig. 11. Because we are interested in accurately measuring colors in the faintest 
regions of the outer disk, we extracted our color profiles from our our uncalibrated B and R 
images, which have slightly superior flatfields to our photometrically-calibrated data; maximum 
amplitude variations of the flatfields are ~1% of the sky in B and ~0.5% of the sky in R. We 
then used comparisons between high signal-to-noise inner disk regions to calibrate the absolute 
colors from the photometrically-calibrated images. The maximum uncertainties expected from 
the combination of sky subtraction and flatficld errors are overplotted as dotted lines on Fig. 11. 
Although signal-to-noise in the outer disk is low, the profiles extracted at all four positions parallel 
to the major axis show regions with B — R <0.60, suggesting these blue outer disk colors are real. 

Finally, we draw attention to the presence of a faint, thicker, but highly flattened disk of 
unresolved stars visible in our color map surrounding the UGC 7321 disk at |r| <2'0. This 
component has B — R and shows little change in color with galactocentric distance (Fig. 12). 
Based upon the observations of Galactic globular clusters (see Seeker 1995), old, metal-poor stellar 
populations are expected to have B — R >0.8, hence the highest z-height stars in UGC 7321 
appear to represent an "old disk" population (see also below). 

3.2.2. Gauging the Effects of Dust 

Bluing as a function of increasing galactocentric distance is a well-known feature of spiral 
galaxy disks (e.g., de Jong 1996). De Jong (1996) has demonstrated through Monte Carlo 
simulations that in face-on spirals, this observed radial bluing generally cannot be accounted for by 
dust. He shows that for realistic models, dust creates color gradients of less than 0.3 magnitudes 
in B — i?, and argues instead that the observed radial color changes in spirals are indicative of star 
formation progressing radially outwards in disks over time, leading to stellar age and metallicity 
gradients. This picture is consistent with a number of semi-analytic galaxy formation models 
where galaxy disks are built "from the inside out" (e.g.. White & Frenk 1991; Mo et al. 1998). 
Nonetheless, the radial color gradient we observe in UGC 7321 {/S.B — ~ 1) is significantly larger 
than typical gradients in the de Jong sample (where typically AS — R ~0.6 magnitudes, even 
when bulge light is included). Moreover, in most spirals observed edge-on, radial color gradients 
are generally found to be small (e.g., Sasaki 1987; Aoki et al. 1991; Wainscoat, Freeman, & Hyland 
1989; Bergvall & Ronnback 1995) or negligible (e.g., Hamabe et al. 1980; Jensen k, Thuan 1982; 
de Grijs 1998). 

Unfortunately, in edge-on spiral galaxies, color gradients become more difficult to interpret 
physically, particularly since internal extinction can be quite significant near the galactic plane. 
The net effect can be both an alteration in the galaxy luminosity profile, and a reddening of the 
observed optical colors which can vary as a function of r and z. Therefore in order to determine 
what fraction of the color gradients we observe in UGC 7321 are due to true age and/or metallicity 
gradients, we must include an assessment of the role of dust in this galaxy. 
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In practice, accurately quantifying the effects of dust on the observed colors and luminosity 
distribution in a galaxy is a complicated problem (e.g., de Jong 1996; Xilouris et al. 1999). 
Moreover, standard assumptions about the nature and distribution of dust in normal giant spirals 
are unlikely to be applicable to low-metallicity LSB galaxies (cf. Han 1992). As a result, the effects 
of internal extinction in a given galaxy are ideally derived using a combination of sophisticated 
radiative transfer models and empirical measurements (e.g., Xilouris et al. 1999). Such modelling 
beyond the scope of the present work, so here we derive some fundamental constraints on the 
effects of dust in UGC 7321 using a very simple model. 

Although some authors have argued that the internal extinction in LSB galaxies may be 
almost negligible (e.g., McGaugh 1994; Tully et al. 1998), these claims have only rarely been tested 
in edge-on systems (e.g., Goad &: Roberts 1981; Kodaira & Yamashita 1996; Bergvall & Ronnback 
1995; Karachentsev 1999). We have already mentioned several lines of evidence that the internal 
extinction in UGC 7321 is low (regularity of the iJ-band isophotes; visibility of background 
galaxies through the disk; low Lfir/Lb ratio; correspondance between the galaxy center in H- 
and S-band). Nonetheless, visible inspection of our WIYN images shows the galaxy is clearly not 
completely devoid of dust, and a number of individual clumps (possibly molecule-rich dark clouds) 
can be seen in our optical images (e.g.. Fig. 5; see also Matthews et al. 1999b). In spite of this, 
the fact that we resolve many of these individual dust clumps instead of seeing a uniform dust 
lane implies that the disk of UGC 7321 is not optically thick. 

Another test of the optical thickness of the UGC 7321 disk comes from a comparison between 
its optical rotation curve (derived from longslit spectroscopic measurements of Ha emission 
from Hll regions), with a rotation curve derived from Hi measurements. Prom optical longslit 
measurements. Goad & Roberts (1981) found UGC 7321 to have a slowly rising, solid-body 
rotation curve throughout much of its stellar disk. If a galaxy disk is optically thick, the rotation 
curve may appear solid body regardless of its intrinsic shape, as an artifact of one's inability to 
observe Hii regions at small galactocentric radii (Goad &; Roberts 1981; Byun 1993). However, a 
rotation curve of UGC 7321 derived from Hi aperture synthesis measurements (Sect. 4.2) confirms 
that the slow rise of the rotation curve throughout the stellar disk is indeed intrinsic (see Fig. 19, 
discussed below). Once again we conclude the disk of UGC 7321 is not optically thick. 

Dust affects the color and luminosity profiles of galaxies through both the scattering and 
absorption of optical photons. However, the maximum color change (i.e., maximum reddening) 
occurs in the pure extinction case. In addition, Byun, Freeman, & Kylafis (1994) have shown that 
scattering effects become decreasingly important with increasing inclination (see also Bianchi, 
Ferrara, & Giovanardi 1996). Since here we are primarily interested in estimating the effects of 
dust on the color gradient, and since we are in the optically thin regime, we can simplify the 
problem by considering a model of extinction due to a foreground dust screen (see Disney, Davies, 
& Phillips 1989). For a given amount of dust, the Screen model produces more reddening than a 
more realistic model of stars mixed with dust, hence it establishes an upper limit to E{B — i?) as 
a function of radius. 
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To evaluate E{B — R), we begin by attempting to attribute as much of our observed color 

gradient as possible to dust. Wc further assume that, to first order, dust extinction in the iJ-band 
is negligible (e.g.. Fig. 2). In the case of zero intrinsic radial color gradient, a dust-free version of 
UGC 7321 would exhibit no radial color gradient in R — H, hence we can use the changes we do 
observe in the R — H color to estimate the amount of extinction at a given radius. We do not 
attempt to reproduce the clumpy nature of the dust in our model (which again, would decrease 
any reddening effects), but assume the dust distribution is roughly exponential in the radial and 
vertical directions (e.g., Wainscoat, Freeman, &: Hyland 1989; Xilouris et al. 1999). In this case 
the extinction along an interval from r to 5r can be expressed as 

SAx,i{r,z) = { ' (4) 



lax 



where ^;^^o is the absorption in magnitudes per unit length at a given wavelength, r = Rma 
corresponds to the edge of the stellar disk, hr^d is the scale length of the dust, and hz,d is the scale 
height of the dust (see also Wainscoat, Hyland, k, Freeman 1989). 

Fig. 13 shows a plot of the "pseudo" R — H color profile of UGC 7321 along the major axis 
interval r =0-120". Here the i?-band data were smoothed to match the resolution of the iJ-band 
observations. Because our H-hand data are not photometrically-calibrated, we have chosen an 
arbitrary normalization such that our "pseudo" R — H color is near r=80". In Fig. 13, we see 
that near r=0, the i2-band light is depressed by roughly 40% relative to the iJ-band, corresponding 
to an i?-band extinction Ar^ RiO.55 magnitudes per kiloparsec. Assuming the extinction law of 
Bouchet et al. (1985), we derive Ab^o=0-89 magnitudes per kiloparsec. 

From Fig. 13 we estimate a scale length of the dust distribution of ~40", which appears 
consistent with the observed dust distribution in our images. Although in giant spirals, it is 
often assumed the dust has a similar scale length to the old stellar disk (e.g., Kylafis & Bahcall 
1987), or even extends beyond it (e.g., Xilouris et al. 1999), this docs not appear to hold in 
UGC 7321. We see an abrupt end to the presence of resolved dust clumps beyond radii |r| 80" 
in both our WIYN images, and in Hubble Space Telescope images obtained by Matthews et al. 
(1999b). Moreover, background galaxies are readily visible through the UGC 7321 disk beyond 
these radii (Matthews et al. 1999b), implying extinction has dropped appreciably. For the vertical 
scale height of the dust, we adopt /i^,d=2"3, which is half the scale height of the old stellar disk 
(Sect. 3.3; see also Paper II). 

Using Equation 4, wc can now derive an extinction-corrected B — R color profile for UGC 7321 
along its major axis (Fig 14a) and offset 1"6 to the north (Fig 146). We see that even the 
extinction-corrected radial color gradient remains pronounced {/S.B — ~ 0.80 mag). This 
demonstrates that a significant fraction of the observed radial color gradient in UGC 7321 appears 
to be due to stellar population and/ or metallicity gradients. 
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3.2.3. Discussion and Interpretation 

The presence of the observed radial disk color gradients in UGC 7321 has several implications. 
Regions with B — R < 0.8 are too blue to be accounted for by old, metal-poor stellar populations 
(e.g., Seeker 1995), but instead must result from a mixture of ages. Without resolving these 
populations, we cannot uniquely constrain the mix of stars present. However, the bluest colors we 
reliably measure in UGC 7321 (B — R < 0.6) are comparable, for example, to young regions in 
the moderate metallicity giant Magellanic galaxy NGC 4449 (Bothun 1986) and to the predicted 
colors of the younger outer regions of a smoothly evolving galaxy with a high halo spin parameter 
and an age of ~10 Gyr, as in the models of Jimenez et al. (1998). In either case, it is likely that 
youthful stars are a significant contributor to the blue colors of the outer disk of UGC 7321 (see 
also Bell et al. 1999). This assertion is strengthened by the presence of Ha emission in the outer 
disk (Fig. 7) and the existence of resolved supergiant stars at these radii (Matthews et al. 1999b). 

We note that edge-on, superthin LSB galaxies like UGC 7321 are especially valuable for 

studies of the outer regions of disks, as they may contain some of the most pristine, unevolved 
galaxian disk environments observable in the nearby universe. However, because these disk 
regions are so optically thin, they would become very difficult to detect in galaxies viewed at 
lower inclinations. For example, if viewed near face-on, the portions of the UGC 7321 disk where 
B — R <0.6 would have an observed 5-band surface brightnesses of only >27.6 mag arcsec"^. Such 
faintncss levels are only rarely reliably achieved in galaxy imaging surveys due to combinations of 
short integration times, scattered light, limited detector ficlds-of-view, and flatfield uncertainties 
(cf. Morrison, Boronson, k. Harding 1994; Lequeux et al. 1996). 

In spite of the potentially very young outer disk regions of UGC 7321, its moderately red 
central disk (where B — R ~1.2) and its somewhat thicker "old disk" (where B — R ~1.1) both 
indicate that UGC 7321 is likely to contain stellar populations with ages in excess of 10 Gyr (see 
Jimenez et al. 1998). As demonstrated above, accounting for a reasonable amount of reddening 
due to dust does not affect this conclusion. The red stellar components in UGC 7321 imply that 
this is not a young galaxy presently undergoing its first epoch of star formation. 

Our observed colors and color gradients in UGC 7321 are similar to those found in other LSB 
spirals by de Blok, van der Hulst, & Bothun (1995; see also Impey & Bothun 1997). A consistent 
model for UGC 7321 is then that of a disk galaxy which has made stars slowly and is now seen to 
be under-evolved relative to typical giant spirals (see also O'Neil et al. 1997; Jimenez et al. 1998). 
This would account for the low gas metallicity of UGC 7321 (Goad & Roberts 1981), the large Hi 
mass fraction (Sect. 4), and the low density of stars (i.e., the low optical surface brightness disk). 
Thus while LSB disks like UGC 7321 may be unevolved galaxies, we concur with the suggestion of 
Jimenez et al. (1998) that LSB galaxies are not necessarily young. 

In order to preserve its relatively strong radial color gradients over many Gyr, it appears that 
viscous evolution in the disk of UGC 7321 has been minimal. In giant spirals, viscous evolution 
(which redistributes angular momentum in the disk) is often argued to be responsible for the 
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relatively smooth exponential luminosity profiles of most disk galaxies (e.g., Lin & Pringle 1987). 

This process is also expected to partially smooth color gradients due to stellar age and metallicity 
changes with radius (e.g., Firmani, Hernandez, & Gallagher 1996). However, such a mechanism 
may not be able to work efficiently in galaxies like UGC 7321, due to both the low disk surface 
density (which considerably extends dynamical timescales) and the slowly rising rotation curve 
[which produces limited shear in the inner disk regions (Lin k. Pringle 1987; see also Matthews 
&: Gallagher 1997)]. The gradients in the H-hand major-axis profile of the UGC 7321 disk 
(Sect. 2.1) and the deviation of the azimuthally-averaged profile from a pure exponential disk at 
small r (Sect. 2.2.3) may be additional signs that viscous evolution has been limited in this galaxy 
compared with normal spirals. 

Finally, we note that at intermediate galactic radii, the eastern side of UGC 7321 has a 
slightly bluer mean color than the western side. This result is insensitive to reasonable sky 
subtraction and flatfield errors. One explanation could be that we are seeing the projection of a 
spiral arm(s), although it remains uncertain how much spiral structure a thin disk like UGC 7321 
might actually sustain, since spiral structure is believed to be an efficient disk heating mechanism 
(e.g., Lacey 1991 and references therein). Moreover, at least some self-gravity is needed for spiral 
arm formation, but if the disk of UGC 7321 were completely self-gravitating, it would be highly 
unstable to "firehose" instabilities and hence to vertical thickening (Zasov, Makarov, &; Mikailova 
1991; Matthews 1998). An alternative is that there are asymmetries in the stellar distribution of 
UGC 7321 due to patchy star-formation (see also Gerritsen & de Blok 1999). 

3.3. Vertical Color Gradients 

The presence of vertical color gradients in galaxies offers an important clue toward the 
dynamical evolution of galaxian disks, since such gradients are predicted to occur as a consequence 
of dynamical heating processes (e.g.. Just, Fuchs, & Wielen 1996). Unfortunately, vertical color 
gradients are very difficult to measure in most edge-on spirals due to dust, contamination from the 
bulge component, and the effects of atmospheric seeing. To date, only a few analyses of edge-on 
spirals have uncovered non-negligible vertical color gradients (e.g., Wainscoat, Freeman, k. Hyland 
1989; Bergvall k, Ronnback 1995), while many workers have reported such gradients to be very 
small or negligible (e.g., Jensen k, Thuan 1982; de Grijs, Peletier, k van der Kruit 1997). 

Fig. l^a-h illustrates the vertical color profiles of the UGC 7321 disk at various galactocentric 
distances. At r = 0, we see very little z color gradient: B — R ^1.2 at a range of z values 
(Fig. 15a), although a slight asymmetry is visible. At r = 6" (Fig 156), we see the addition of a 
peak with B — R ~1.45. Note its displacement from the disk midplane; this peak is in the vicinity 
of the red nuclear feature discussed in Sect. 3.1. At r=10" (Fig 15c), the vertical color gradient is 
again fairly flat {B — R ~1.2), with a very slight bluing visible near z=l".^. At r =0'5 (Fig 15(i), 
the high z regions of the disk have nearly the same color as those at r=0 (this is due to the "old 
disk" of red stars discussed above), but at small z values, the disk has become much bluer than 
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at smaller galactocentric radii. Finally, with further increasing radius, while the color at high 
z continues to stay nearly constant, the color at small z grows increasingly bluer (Fig. 15e-/i), 
consistent with the radial bluing observed along the major axis of the galaxy in Fig. 11. 

Vertical color gradients of the type we observe over much of the disk of UGC 7321 (i.e. redder 
with increasing z height) are predicted to occur due to dynamical disk evolution in which older 
stellar populations acquire higher velocity dispersions over time due to heating processes (e.g., 
Fuchs h Wielen 1987; Just, Fuchs, k, Wielen 1996). UGC 7321 represents one of the few examples 
of this type of vertical color structure being directly observed in a galaxy disk. Thus even the 
most dynamically cold examples of nearby galaxy disks appear to have undergone some dynamical 
heating. Nonetheless, the exact mechanism by which older stars are dynamically heated and how 
they are redistributed as a function time remain uncertain (see review by Lacey 1991) and different 
mechanisms may operate in different galaxies, depending upon factors like environment, degree of 
self-gravity, and the size and number of molecular clouds. The relative simplicity of UGC 7321, 
and its status as a relatively unevolved galaxy make this an ideal system to help place important 
constraints on these issues. We explore this further in Paper II. 

Using the extinction parameters derived in Sect. 3.2.2, we can also estimate the effects of 
dust reddening on our observed vertical color gradients. We adopt for a scale height of the dust 
half the scale height of the old stellar disk (Xilouris tt al. 1999). Using the stellar scale height 
of the old disk derived in Paper II, this yields hz^d='^"-3- Fig. 16 shows an example of one of our 
vertical color profiles corrected for dust reddening. We see that at the galactocentric radii where 
the vertical color gradients become the strongest, correction for dust reddening acts to slightly 
increase A{B — R) along the vertical direction. 

4. HI Observations 

4.1. HI Pencil Beam Mapping 

In June 1997 we used the Nancay Decimetric Radio Telescope to obtain a 7-point pencil-beam 
map of UGC 7321 in the 21-cm line of neutral hydrogen (H i). At the declination of UGC 7321, the 
Nangay Radio Telescope has a FWHM beam size of approximately 4'E-Wx23'N-S. Other details 
regarding the Nangay telescope may be found in, e.g., Matthews, van Driel, k, Gallagher (1998). 

UGC 7321 was observed at 7 different telescope pointings, including a position corresponding 
to the optical center of the galaxy, and positions offset 2', 4', and 6' east and west of center, 
respectively. The observations were obtained in total power (position-switching) mode using 
consecutive pairs of two-minute on- and two-minute off-source integrations. Total integration 
times were 2-3 hours at each pointing. For these observations, the auto correlator was divided into 
two pairs of cross-polarized receiver banks, each with 512 channels and a 6.4 MHz bandpass. This 
yielded a channel spacing of 2.6 km s~^, for an effective velocity resolution of ~3.3 km s~^. 
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Fortuitously, UGC 7321 is oriented along nearly an E-W line, so the Nangay telescope 
provides sufficient spatial resolution along this direction to obtain crude information about the 
Hi distribution and kinematics of the galaxy. The individual spectra obtained at each telescope 
pointing are shown in Fig. 17 and the resulting global spectrum in Fig. 18. Note the global profile 
appears quite symmetric. Table 1 gives our derived global H i parameters. Errors were computed 
following Matthews, van Driel, k. Gallagher (1998). Our values are in good agreement with the 
recent singlc-dish Hi observations of UGC 7321 obtained by Haynes et al. (1998) using the Green 
Bank 140-ft telescope. 

Our Nangay observations show that the Hi gas in UGC 7321 clearly extends beyond the 
optical galaxy, confirming the results of Hewitt, Haynes, & Giovanelli (1983). Significant flux 
is detected at the pointings both 4'E and 4'W of the galaxy, giving a lower limit for the for 
H i-to-optical diameter ratio of Djji/Dopt >l-25. In addition, we detect a small amount of Hi flux 
(~1.1 Jy km s~^) in our observation 6'W of the galaxy center. From the Nangay data alone, it is 
difficult to assess whether this flux is due to a sidelobe contamination or real extended emission, 
since the strength of the Nangay telescope sidelobes varies significantly with the hour angle of the 
source (Guibert 1973). 

4.2. HI Aperture Synthesis Data 

After our Nangay observations were obtained, we had the opportunity to examine H I aperture 
synthesis data of UGC 7321 obtained by Rots, Roberts, & Goad (private communication) using 
the Very Large Array (VLA). These observations were obtained in C array in 1981 using 18 
antennas, hence they do not achieve the same sensitivity limits as modern VLA data, but they still 
offer a useful complement to our Nangay dataset. Thirty-one independent velocity channels were 
used for the observations, with a velocity resolution of ~10.3 km per channel. The FWHM 
beamwidth was ~12". 

Using the VLA data, we derived a position- velocity diagram for UGC 7321. By assuming 
that at each point along the major axis the maximum velocity traces the rotation, we have derived 
the rotation curve shown in Fig. 19. This figure confirms the slowly-rising, solid body nature 
of the rotation of UGC 7321 throughout its stellar disk, as was first seen in the optical rotation 
curve of Goad Sz. Roberts (1981). Although beam-smearing may be expected to slightly decrease 
the amplitude of our derived H I rotational velocities, this effect is not expected to exceed a few 
kilometers per second, or to be significant beyond 2-2.5 beam diameters from the disk center 
(e.g., Swaters 1999), hence it cannot explain the shallow shape of the inferred rotation curve. 
Finally, wc note that Fig. 19 suggests that the rotation curve of UGC 7321 does not begin to 
flatten until near the edge of the stellar disk. This is consistent with UGC 7321 being a galaxy 
whose dynamics are dominated by a dark halo even at small galactocentric radii (e.g., Matthews 
1998). However, deeper, more sensitive Hi aperture synthesis measurements would be valuable for 
further constraining the detailed shape of the outer rotation curve. 
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The integrated Hi map of UGC 7321 derived by Rots, Roberts, & Goad (private 

communication) confirms a slight extension of the H I on the west side of the galaxy compared 
with the cast side, consistent with our Nancay data. It also permitted us to measure an Hi 
diameter of Dhi ~7'1 for UGC 7321 at the limiting Hi column density of the observations 
{Nhi ~ 7 X 10^0 atoms cm ). Thus Dni/Dopt >l-2, in agreement with our Nangay lower limit. 



4.3. Discussion 

It is somewhat difficult to accurately compare the Hi extent of UGC 7321 with other 
spirals since the commonly-quoted DH1/D25 ratio (where D25 is measured at a face-on-corrected 
isophote) in not particularly meaningful for a galaxy like UGC 7321, whose face-on central surface 
brightness is only ~ 23.4 mag arcsec"^. Moreover, we cannot unambiguously translate our limiting 
observed Hi column density into an Hi surface density due to projection effects, thus it is difficult 
to accurately measure Dhi at some canonical surface density (e.g., 1 pc~^). 

In terms of optical scale length, the Hi in UGC 7321 extends to at least 5.5hr. Although 
this value may increase with more sensitive VLA observations, such an extent is still quite normal 
for both high and low surface brightness late-type spirals having similar rotational velocities to 
UGC 7321 (see Fig. 8 of de Blok, McGaugh, & van der Hulst 1996). 

If we assume both the stellar and the Hi disks of UGC 7321 are optically thin, and hence 
both similarly enhanced in surface brightness due to their edge-on projection, then we see that the 
maximum Hi extent of the disk relative to the maximum optical extent is typical for Sd spirals 
(e.g., Hewitt, Haynes, & Giovanelli 1983). We can more directly compare with two other Sd 
superthins previously mapped in Hi: Dni/Dopt ~1.2 was also found by Cox et al. (1996) for the 
superthin spiral UGC 7170 and by Abe et al. (1999) for the superthin IC 5249, both at similar 
limiting Hi column densities and optical surface brightnesses. 

Ignoring internal extinction effects, the ratios of UGC 7321, UGC 7170, and IC5249 are 
2.6, 3.4, and 2.3 (in solar units), respectively, and the superthin spiral ESQ 146-014 studied by 
Bergvall & Ronnback (1995) has =^~2.8. If we assume all have internal extinctions similar to 
UGC 7321 {AB,i ~ 0.89 mags; Sect. 3.2.2), the implied ^^jj^ ratios are still significantly larger than 
typical values for normal Sd spirals derived by Roberts &; Haynes (1994; =^^~0.63), although 
not as extreme as some of the unevolved dwarf galaxies studied by van Zee et al. (1995) or some 
of the extreme late-type LSB spirals in the sample of Matthews & Gallagher (1997) which have 
•^^^ as high as 10 (see also Salzer et al. 1991; Matthews, van Driel, & Gallagher 1998). 

Van der Hulst et al. (1993) suggested that low gas surface densities are responsible for the 
low star-formation efficiencies in LSB galaxies (see also de Blok, McGaugh, & van der Hulst 1996; 
Gerritsen & de Blok 1999). Although we cannot unambiguously recover the radial Hi surface 
density distribution in UGC 7321 due to its edge-on geometry (cf. Oiling 1996), we can compute 
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a mean Hi surface density within the stellar disk. Adopting the definition T,hi =M.Hi/'^Rlpt: 
where Ropt is the linear optical diameter of the disk (e.g., Roberts & Haynes 1994) we find 
^H/=5-3 Mq pc~^. This value is considerably lower than the median of this quantity found by 
Roberts & Haynes (1994) for Scd/Sd spirals in the UGC catalogue: I]ji-/=9.80 M.q pc~^, but is 
consistent with the T^hi ratios of several of the LSB galaxies in the sample of de Blok, McGaugh, 
& van der Hulst (1996). 



5. Summary 

We have presented B- and i?-band imaging and photometry, near-infrared iJ-band imaging, 
narrow-band Ha imaging, and H I 21-cm line measurements of the nearby, edge-on Sd spiral galaxy 
UGC 7321. 

UGC 7321 is a "superthin" galaxy, with an extremely dynamically cold stellar disk and 
no discernible bulge component, even in the near-infrared. In spite of its edge-on orientation, 
UGC 7321 is visibly quite diffuse and it is clear that its intrinsic optical surface brightness is quite 
low. The dust content of UGC 7321 also appears to be small, hence we argue that the "anemic" 
appearance of the galaxy results from the low-level of current star formation rather than severe 
internal extinction. 

UGC 7321 exhibits significant B — R color gradients in the radial direction: measured from the 
disk center to its edges edge, A{B — R) ~1.05 magnitudes along the galaxy major axis. Dust alone 
cannot explain the large gradient. Using a simple extinction model we find A{B — R)cor ~0.80 
after correction for internal reddening. This is somewhat larger than the color gradients typically 
observed in normal giant spirals (cf. de Jong 1996), and suggests significant stellar population 
gradients in the disk of UGC 7321. 

The outermost disk regions of UGC 7321 have B — R <0.6, suggesting they are composed 
of stellar populations that include a significant fraction of young stars. However, UGC 7321 also 
contains a population of old stars with B — R >1.1, indicating it is not a young or recently- formed 
galaxy. The rather strong radial segregation of these populations suggests that the galaxy has 
evolved quite slowly and perhaps that viscous evolution has not operated efficiently in this system. 

UGC 7321 also exhibits appreciable vertical color gradients: AB — i? as large as 0.45 
magnitudes was measured parallel to the minor axis. This is a reflection of a concentration of 

young blue stars along the galaxy midplane, and a population of the older, red stars at larger scale 
heights. This type of age segregation is predicted to occur due to dynamical heating processes 
in spiral disks, but UGC 7321 represents one of the few examples of it being directly observed in 
an external galaxy. This implies that even dynamically cold disks like UGC 7321 have undergone 
some dynamical heating. 



The stellar light distribution of the UGC 7321 disk cannot be characterized by a single 
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exponential function. A light excess over the prediction of an exponential model is seen at small 
radii, and in both optical and NIR wavelengths, the disk appears to contain distinct surface 
brightness "zones" between which the slope of the brightness profile changes. This may be an 
additional signature of minimal viscous evolution. It also brings into question the prediction of 
some semi-analytic galaxy models (e.g., Dalcanton, Spergel, k. Summers 1997) that exponential 
stellar disks are a natural product of the disk formation process. 

In spite of its regular, organized disk, many of the global properties of UGC 7321, including 
its luminosity. Hi content, rotational velocity, and "^^^ ratio (Table 1), as well as its slowly-rising, 
solid-body rotation curve (Fig. 19) are more typical of a dwarf irregular galaxy than a normal 
Sd spiral. The origin of such vastly different disk morphologies in an otherwise similar physical 
parameter space remains unclear, but it may place an important constraint on galaxy formation 
and evolution models. 

Together the properties of UGC 7321 suggest that it is an under-evolved galaxy in both a 
dynamical and in a star-formation sense. Nonetheless, this galaxy clearly demonstrates that even 
seemingly simple "pure disk" galaxies like the superthins are highly complex systems. 

We are grateful Arnold Rots and Mort Roberts for providing us with the fully-calibrated 
VLA Hi data of UGC 7321. We also thank Mike Merrill for assistance with the IRIM observations 
at Kitt Peak and Wanda Ashman for creating the artwork for the Appendix. This research 
was partially funded by the Wide Field and Planetary Camera 2 (WFPC2) Investigation 
Definition Team, which is supported at the Jet Propulsion Laboratory (JPL) via the National 
Aeronautics and Space Administration (NASA) under contract No. NAS7-1260. The Nangay 
Radio Observatory is the department Unite Scientifique Nangay of the Observatoire de Paris and 
is associated with the French Centre National de Recherche Scientifique (CNRS) as the Unite de 
Service et de Recherche (USR), No. B704. The Nangay Observatory also gratefully acknowledges 
the financial support of the Region Centre in France. This research made use of the NASA/IPAC 
Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aeronautics and Space Administration. 

A. Estimating the Effects of Inclination on the Observed Central Surface 

Brightness of a Disk 

When deriving central surface brightness measurements of galaxies, it is desirable to convert 
measured values to face-on values in order that galaxies observed at different inclinations 
can be meaningfully compared. However, for highly-inclined galaxies, the approximation that 

IJ' face-on = ^^observed + 2.51og(a/6) (e.g.. Freeman 1970) breaks down, and alternative corrections 
must be used. Here we derive a simple formula for this purpose, using a "thin disk" approximation. 

In our model we assume the galaxy disk is optically thin and has an exponential brightness 
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profile along the r direction of the form L(r) = Loe~'^/^^, where Lq is the luminosity volume 
density and hr is the disk scale length. We also assume the disk has a finite but non-negligible 
thickness of 2hz, and a constant brightness along the z direction (see Fig. 20). For such a disk, the 
central surface brightness if viewed edge-on is expressed as 

L(0),=90 = Lo / e-'-Z'^'-dr = 2Lo/i.; (Al) 

if viewed face-on, the central surface brightness would be simply L(0)i=o = 2Lohz- 

Now suppose our model disk is viewed at some arbitrary angle i' < 90° through a line of sight 
along a vector x (see Fig. 21). We define i as the angle between the line of sight and a normal to 
the disk. Now the integrated brightness of the disk along a path from a; = to x = x^ax becomes 

L{0)i=i, = 2Lq / e^^'^'/'^'-sim dx. (A2) 
Jo 

Taking the integral, and making the substitution 

Xmax — hz/cosi, (A3) 

we find that 

L{0)i = 2Lohr[l - e"'*-*^'^'/'*-]. (A4) 

To apply this approximation to a real galaxy, one can substitute the measured radial and vertical 
scale heights for hr and hz respectively, provided that the conditions hz « hr and 

(90 - i') > tan-l {hz/Xrr^ax) (A5) 



both hold. 
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Fig. 1.— H-hand image of UGC 7321 obtained with IRIM on the KPNO 2.1-m telescope. Field 

size is ~4{6x l^O. East is on the right of the image, north on top. A number of foreground stars and 
cosmic rays have been cleaned from the image. Note the complete absence of a bulge component, 
even in the near-infrared. 

Fig. 2. — Contour map produced from the i?-band image of UGC 7321. The image has been 
stretched in the z direction for illustration purposes. The contours are spaced in units of 20 counts. 
Note the regularity of the contours, even near the galaxy center, suggesting dust absorption is 
minimal in the iJ-band. 

Fig. 3. — //-band major axis profile of UGC 7321, averaged over a 3-pixel-wide strip. The axes 
are distance along the galaxy major axis, in arcseconds, versus the logarithm of the intensity, in 
arbitrary units. The arrows mark the locations of changes in slope of the brightness profile. 

Fig. 4. — ii-band image of UGC 7321 obtained with the WIYN telescope. Image size is ~6'5xl'0. 
East is on the right, north on top. Seeing was ~0''6 FWHM. 

Fig. 5.— Close-up of a 1(8 long portion of the i?-band WIYN image of UGC 7321. The optical 
center of the galaxy corresponds to the large black patch on the righthand side of the image. White 
areas show dust clumps in the disk. Note the absence of a true dust lane. 

Fig. 6. — Close-up of the western edge of the UGC 7321 disk in the R band, showing a background 
galaxy visible though the disk, at approximately 2'.5 west of the disk center. 

Fig. 7. — Continuum-subtracted Ha image of UGC 7321 obtained with the WIYN telescope. Scale 
and orientation are as in Fig. 4. 

Fig. 8. — Detail of the Ha emission structure centered at ~ 8" east of the galaxy center. This 
image section is roughly 50" long. Note the gap in the emission at the center of this frame and the 
faint plumes extending out of the disk midplane. 

Fig. 9. — Azimuthally-averaged B-band surface brightness profile of UGC 7321 obtained from 
fitting a series of 9 concentric ellipses. Axes are distance from the galaxy center in arcseconds, 
versus average S-band surface brightness in magnitudes per square arcsecond. The dashed line 
represents the best exponential fit to the data. No inclination or extinction corrections have been 
applied. 

Fig. 10. — B — R color map of the inner 3'8 of UGC 7321 obtained from our WIYN data. Seeing 
in both the B and R data was ~0''6. Black (near the nuclear region) corresponds to _B — i? ~ 1.5; 
dark grey (seen surrounding the nuclear region): B — R ~1.2; mottled black and white (seen in the 
thicker disk surrounding the galaxy): B — R Ril.l; white (as seen along the disk major axis and in 
the outer disk regions): 5 — < 0.85. 

Fig. 11. — (a) B — R color along the major axis of UGC 7321, averaged over a 12-pixel-wide strip. 
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The data have been smoothed by a factor of 15 along the r-direction. The axes are distance from 

the galaxy center in arcseconds versus B — R color in magnitudes, corrected for Galactic foreground 
extinction, (b) Same as (a), but offset 1"6 north of the major axis, and running through the very 
red nuclear feature (see Sect. 3.1). (c) Same as (a), but offset 1".2 south from the major axis, (d) 
Same as (a), but offset 4''6 south from the major axis. The very red feature near r = —135" is due 
to a background galaxy seen through the disk. 

Fig. 12. — Same as Fig. 11(a), except the profile was extracted 10" north of the major axis. This 
cut runs through the thicker "old disk" of stars visible in Fig. 10. Note there is very little color 
gradient as a function of r. Beyond r = ±75" our extracted profile no longer traces the "old disk" 
due to the curvature of the galaxy isophotes. 

Fig. 13.— "Pseudo" R-H color along the major axis of UGC 7321, in the interval r < 120". The 
R data were smoothed to the iJ-band resolution, and the R and H intensities have been arbitrarily 
normalized such that R — H near r=80". 

Fig. 14. — (a) Same as Fig. 11(a), but with internal reddening correction applied. (6) Same as 
Fig. 11(6), but with internal reddening correction applied. 

Fig. 15. — B — R vertical color profiles extracted parallel to the minor axis of UGC 7321. The 
profiles are averaged over 15-pixel- wide strips and smoothed along the z direction. Axes are distance 
from the galaxy midplane, in arcseconds, versus B — R color, corrected for Galactic foreground 
extinction, (a) Minor axis profile (r=0); {h) r=5" east; (c) r=10" east; {d) r=0{5 east; (e) r=l'0 
east; (/) r=l{25 east; {g) r=l'5 east; {h) r=2{0 east. 

Fig. 16. — Vertical B — R color profile at r=l{25 east of the disk center, corrected for internal 
reddening. The uncorrected profile is overplotted as a dashed line. 

Fig. 17. — Hi pencil beam map along the major axis of UGC 7321 obtained with the Nangay Radio 
Telescope. The center panel shows the spectrum obtained at the optical center of the galaxy, and 
the remaining panels show spectra obtained at pointings offset west and east of center, respectively, 
at multiples of one-half beamwidth. Axes are heliocentric radial velocity, in km s~^xlO~^, versus 
flux density in millijanskys. 

Fig. 18. — The global Hi profile of UGC 7321. Data from the seven telescope pointings shown in 
Fig. 17 were combined following Matthews, van Driel, & Gallagher (1998). Axes are heliocentric 
radial velocity, in km s~^, versus flux density in millijanskys. 

Fig. 19. — Rotation curve of UGC 7321 derived from the VLA Hi aperture synthesis data of 
Rots, Roberts, & Goad (unpublished). Axes are distance along the major axis, in arcseconds, 
versus heliocentric radial velocity, in kilometers per second. The horizontal dashed line denotes the 
systemic velocity measured from our Nangay Hi pencil beam measurements. The arrows indicate 
the edge of the stellar disk in our WIYN i?-band observations, and the crosses denote the limiting 
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extent of the optical rotation curve derived by Goad &; Roberts (1981). The close agreement between 
the location of Vmax of the rotation curve and the global Hi profile horns in Fig. 18 suggests are 
radial velocities are accurate to ~10 km s~^. 

Fig. 20. — Schematic drawing illustrating a thin disk inclined at some angle i to the line of sight. 
Fig. 21. — Blow-up of the central portion of the disk shown in Fig. 20. 
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Table 1. UGC 7321 Parameters 



General: 

a : 12''15™03^ (1950.0) 
5 : +22°49'03" (1950.0) 
Hubble type: Sd IV 
Distance = 10 Mpc " 
Ab = 0.04 magnitudes * 

Measured Optical Parameters: 

niB = 13.84±0.03 
niR = 12.85±0.03 
Dopt = 5^6 
PA = 82° 
a/h = 10.3 " 

hr = 44" ± 2" = 2.1 kpc f 
//^(O) = 21.6 mag arcsec"^ ^ 
A«B(0)ej; = 22.0 mag arcsec"^ 
jjg = 23.9 mag arcsec"^ ' 
B-R = 0.99 



Table 1. (continued) 



Derived Optical Parameters:-' 

Mb = -17.05 
Mr = -17.70 

i.B = l. 0x10^1.0 

i = 88° 

Aopt = 16.3 kpc 
yU5_j (0) = 23.6 mag arcsec"^ ' 
jjg^i = 27.6 mag arcsec"^ "* 
Lfir/Lb=OM " 

Measured HI parameters: 

14^20 = 233±7 km s'^ 

14^50 = 216±4 km s"! 

Vh = 407.2±2.2 km s"! 

/ Sdv = 47.6±4.0 Jy km s"! 

Dhi=T.I (at Nhi ~ 7 X 10^° atoms cm'^) 

Derived HI parameters: 

Mhi = I.IxIO^Mq 
■Mhi/Lb = 1.1 (solar units) 
Dni/Dopt >1.2 



^From Mattliews et al. (1999b) 
'='From Burstein & Heiles (1984) 

■^Observed magnitudes and luminosities have been corrected for Galactic foreground extinction only; 
isophotal diameters and surface brightnesses are measurec? values-i.e., uncorrected for inclination. 

■^Maximum observed angular extent of the galaxy, corresponding to the observed 25.5 mag arcsec"^ 
5-band isophote, uncorrected for inclination. 

'^Disk axial ratio measured at the observed 25.0 mag arcsec"^ i?-band isophote. 

^Disk scale length measured from the 5-band image. 

^Measured central surface brightness from 5-band image, uncorrected for inclination and internal extinc- 
tion. 

'^5-band central surface brightness obtained from extrapolation of an exponential fit to the outer disk, 

uncorrected for inclination and internal extinction. 

'Mean 5-band surface brightness over the observed extent of the disk, uncorrected for inclination. 

J Derived optical quantities have been corrected for Galactic and internal extinction, as derived in 
Sect. 3.2.2. ^ 

'^Linear diameter at limiting observed 5-band isophote of 25.5 mag arcsec"^. 
'Predicted face-on central surface brightness, based upon the measured quantity //^(O). 
""Predicted mean face-on surface brightness. 

"Far-infrared luminosity was derived from co-added IRAS 100//m fiux from Sage (1993) and IRAS 60//m 
fiux from the NED database 
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